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The propagation characteristics of uniform, cylindrical plasma columns of circular cross 
section in axial, constant magnetic fields are determined. The ratios (plasma wavelength)/ 
(free-space wavelength) and (power flowing in the plasma) /(power flowing outside) are 
evaluated and discussed as functions of the diameter/wavelength ratio and of the plasma 
permittivity, for propagation of circularly symmetrical modes. Brillouin diagrams are also 
derived and compared with the results predicted by the quasi-static approximation. The 
analysis is for a nonattenuating plasma, but formulas giving the attenuation due to electron- 
molecule collisions are also included. 

1. Introduction 

A magneto-ionic cylindrical plasma column in a static axial magnetic field may efficiently 
propagate an electromagnetic wave. This fact may be of iinportancc in interpreting phenom- 
ena like plasma oscillations and atmospherics, and in designing new iiietliods for plasma diag- 
nostics and milhmeter wave generation. 

Whereas varions anthors liave derived tlie characteristic equation both for a straight 
column of circular cross section [Akliiezer et al., 1958] and for an infinite slab [Schumann, 
1958] of uniform plasma, no detailed discussion of the propagation characteristics is given, 
except under the so-called '^quasi-static'^ approximation [Smullin and Chorney, 1958; Trivel- 
piece and Gould, 1959]. 

In such an approximation the propagation wavelength must be much less than the free- 
space wavelength, a condition usuall\' fulfilled only when plasma transverse dimensions are 
much smaller than the free-space wavelength. 

Our analysis, instead, is concerned witli the solution of the exact characteristic equation, 
written for the propagation of circularly synnuetrical modes in a uniform, nonattenuating 
plasma of circular cross section; the permittivity tensor is chosen as constant parameter for 
each solution. We have found it convenient to plot our results as the ratios (plasma wave- 
length) /(free space wavelength) and (power flowing in the plasma) /(power flowing outside) 
versus the diameter/wavelength ratio. 

Using these results we have also derived the plots of angular frequency versus wave number 
for each mode (Brillouin diagrams), assuming constant values for the plasma and cyclotron 
frequencies. These diagrams are snnilar to those obtained using the quasi-static approxima- 
tion, except for a few features, which will be described. 

A final section is devoted to the attenuation formulas, when the electron-molecule col- 
lisions are taken into consideration. 

2. Plasma Dielectric Constant 

The dielectric constant of a uniform plasma in a constant magnetic field £ is a tensor 
Hell, defined by the relationship: 

D, = e,^E^ (1) 

1 Sponsored by the Air Force Cambridge Research Center of the Air Research and Development Command, United States Air Force, through 
its European Office. 
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(2) 



where D is the electric displacement and E is the electric field. 
-> 
If the field B is directed along the positive z axis {i, l = ?>), the en components due to the 

electron motion alone are given, for zero agitation and drift velocities, by the expressions 

[Allis, 1956]: 

r" 2 2 — 1 ^ 

_ _ , _;"\- ij '^ I "« 

eu-..2-a„(e, J., )-.o [_1 , 2,;c.k+i(co^a,„)l"^2M.,+.;(co+a„)|J 

Here z^c is the electron collision frequency (assumed to be velocity independent), co is the angular 
frequency of the applied field (well above the ion gyrofreque ncies, so that the ion contributions 
may be neglected), w^ is the angular plasma frequency {=-^ne}lmeQ, n being the electron den- 
sity, m and e the electron mass and charge, and eo the free space permittivity) and co& is the 
angular cyclotron frequency (=eBlm). 

When z^c<<w^ ^^iid ^c<<(<^— ^&)^ 1^^ above fornnilas become: 

2 2 / -2 \ •>\ "S 
€] — j€, =1H 7, 7^~JVc —^ ^, 

co^ — CO ^{<^b' — ^ ) 



CO CO^ — CO (CO^ — CO j 



• fr 1 -^p • ^p 






(3) 



Since the collision frequency Vc is proportional to the gas pressure, the above inequalities are 
fulfilled, whenever the pressure is sufficient!}^ low; we assume throughout the paper that this 
is always the case. 




Figure 1. Regions of permitted and forhidden ei and 
€3 values. 
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The propagation charactc^ristics will bo first investigated assuming Vc=0 (and then 
^'i =^2 =^^' = 0), a condition which provides propagation without losses. Later we shall 
compute the losses, removing this assumption, but still considering the e^' terms small compared 
to the e ones. 

From (3) one can derive the following relationship between ei, e-z, e^: 

e|^(e,-l)(ei-63). (4) 

It is convenient to assume €i and €3 as the independent variables of our problem. These dielec- 
tric constant components are functions of co, w^ and o^j,; the range of values that ei and €3 may 
have is shown in figure 1. In this figure we have divided the plane of possible values in five 
regions: a separate discussion of the propagation characteristics will be presented for each 
region. In the figure a new convenient frequency ccc, defined as V^^ + <^&; has also been 
introduced. 

3. Propagation in a Cylindrical Column 

We assume a cylindrical plasma of circular cross section in an infinite free space. The 
plasma is uniform and in steady state conditions ; its axis is parallel to the static magnetic field 
direction (s=axis). Our analysis is here limited to the propagation of circularly symmetrical 
modes, namely to the propagation of fields without azinnitlial dependence. 

In a cylindrical orthogonal s^^stem of coordinates (p,(p,z) these fields are [Epstein, 1956] 
inside the plasma: 

E, =AMx,pIR) +A,Mi,p/Ii) ^ 

K=-J ^^ {^ [a f; Mx,p/R)+A, I Mf.p/R)] 

E,J^ [a, ^' Mr,p/R)+A, f^ Mx,p!R)] 

H -_?^ « r^, £.. J^(,^pfJi)+A, f^ J,(.r,p/7?)1 

Aq /Oo L '^1 -^2 J 

.f r^, y J:(.r,p/i?) +A, I J.(^,p/A>)1 



A( 



(5) 



outside the plasma: 



E,=BKo{x,pIB) 



E,= -j~a-Bl-K,{x,plR) 

Af) ^0 



Efp==- 



Ao i/Q 



H,=j^CK,(x,plB) 
IL=~^C^K,{x,pIB) 



Xq Zq Xq 



II.=-j 



.2tB 1 



Xq 



-^B-K,{xopIR) 

Zjq Xq J 



(6) 
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where 



i^=;r=radius of the plasma cyhnder 



a=7^=free space Xo to guided X^ wavelength ratio 



Zo==free space characteristic impedance 



e7o,</b-K'o,S'i=Bessel functions, according to McLachlan symbolism [McLachlan, 1955]. 



Ao 



Xq- 



2tR 



{a^- 



■1)5- 



^1,2^^^ 



€3(€i — Q:^) € 



9 1,2 = 



6i(a^-l) 
e,dl2(a'-l)-e,{e-a') 



g,[_,i^,+,M^^] 



ae2 



(7) 



It is desired that the fields outside the plasma decay radially to zero at infinity. Since Kq 
is analogous to a negative exponential only when the argument is real positive, then Xq must be 
real, and consequently a^>l. 

The condition that the tangential fields be continuous at the plasma boundary gives the 
following transcendental equation, which determines the unknown quantity a^: 



where 



F(Xo,82,^3)=-rF(xo,d„es), 
1 JiidXo) , Ki{Xo) 



F(Xo,8,e^) = 



8 Jo{8xo) Koixp) ^ 
es Ji{8xo) K^JXqY 
8 Jo{8XoyKo{xo) 



(8) 



(9) 



■^2 l + Vl+4e3(€i-l)aV(6-63)(a^-l)^ 



(10) 



Equation (8) is the characteristic equation of our propagation problem. Its solutions will 
be discussed in the next section. 

4. Wavelengths in Plasma Guided Propagation 

For the solution of (8) and a survey of the results it is important to know the behavior 
of i^(a;o,5,e3) as a function of Xq, for different 8^ values. This is plotted in figure 2 and in figure 
3 for positive and negative €3 values respectively. Important points on the Xq axis are X^^\ 
fi^^\8), and v^^\8,e^), which are the mth nonzero solutions of the equations: 

J„(X«)=0 

1 Ji(m<'"'5)_ i?i(M""') 






€3 Jiiv 






8 JoC/^'S) Ko(/"'>)J 



(11) 
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If 8'-^+<^, then 



lid 



>xr/8-^o. If 52^0, 



and v^'^\es>0)->Xr/d 



i/(^)(e3<0)-^Xi"' 'V5-^°° vvhen rn^l, whereas v^'\es<CO) remains finite. 

It is convenient to represent the solutions of (8) as relative wavelengths Xg/Xo=l/a: versus 
(i/Xo for different sets of €1 and €3 values. For this purpose we have proceeded as follows: 

1. From a set of real ei, €3 and a^ values the corresponding dl, 8l and r values are computed. 

2. The Xq solutions of (8) are determined. More than one Xq value may be found; in this 
case for each solution there is a different field configuration, that is called propagation mode. 
Since F(xo,-^8j es)=F{xo, — 8, e^), no sign ambiguity results in taking the square root of 8^. 

3. The corresponding d/Xo are obtained from the equation: 



d 



Xq 



^0 T'^a^—l 



(12) 
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Figure 2. Behavior of ¥ versus \o for e3>0 and 
different 6- values. 
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Figure 3. Behavior of F versus Xq for e3<^0 and 
different 8^ values. 



The computations are repeated for the sauie €1 and e-^, but for different a" values between unity 
and infinity. 

The analysis of (8) for arbitrary ei and €3 values would be unnecessarily general. We have 
then restricted our investigation onl\^ to the sets of values, which may be found in the five 
significant regions shown in figure 1. The following results have been obtained. 

1. In region I (see appendix A and fig. 4) no propagation of the mth mode is possible for 
d/\o values smaller than: 



/dv^'^x^r 

\><0/c TT Le3|€i — 1 



^1 



1 — lH-Ve3(€l— l)(€i— €3) 



]* 



(13) 



At this cutoff condition Xg = Xo; for larger d/Xo values propagation takes place with a \g con- 
tinuously decreasing towards the infinite plasma value: 



«o [6i + V(«i-l)(ei-63)]^ Vei + ea 



(14) 



This limit represents the well-known propagation of a TEM plane uniform wave, with a right- 
hand circular polarization. In fact, at this limit, the Ep and E^ components become equal 
and 90° out of phase. These curves have the same shape, as those for propagation in the 
dielectric waveguides; moreover, by setting 6i = e3=e into (13) and (14) we obtain exactly the 
cutoff condition and the free space propagation wavelength for dielectric waveguides, TEom or 
TAIom modes. 
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Figure 4. Wavelength ratios Xg/Xo versus d/Xo m 
region I. 






, 




' 




€r--2 


^3 = 0.5 






\ 






> 


AJ 






/ 


I SOL y 


' 


n soL^_/ 


/ 



0.6 0.6 



Figure 5. Wavelength ratios Xg/Xo versus d/Xo m 
region III. 
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Figure 6. Wavelength ratios Xg/Xo versus d/Xo in 
regions I V and V. 



2. Ill region II no propagation is ever possible (see appendix A). 

3. In region III (see appendix A and fig. 5) an infinite number of propagation modes exist. 
All the Xg/Xo versus d/\o curves start from the origin, the initial slopes being 7r/v^^\8^), where 
5i = — €3/61. Unity wavelength ratios are attained at the same (d/\o), defined by (13). How- 
ever, when the \g curves have the hook shape shown in figm^e 5, propagation takes place also 
at larger d/\o ratios. In this case we have two wavelength ratios for each d/\o condition; in 
general the group velocities of these two possible waves are opposite. The group velocity is 
zero at the maximum d/\o (points A). 

4. In region IV (see appendixes B and C and fig. 6), as in region III, an infinite number of 
propagation modes exist and all the X^/Xq versus d/\o curves start from the origin, the initial 
slopes being Tr/v^'^^d^). Here, however, propagation is possible over the entire cZ/Xq range; 
the curves have positive slopes and approach, as d/\o-^<^ , an horizontal asymptote which is the 
same for all modes, except for the mode of the lowest Xq solution. 

5. In region V (see again appendixes B and C and fig. 6) we have at most one X^/Xq versus 
d/Xo curve for each pair of €1,63 values. The curve has the same behavior as the curves found 
in region IV. The initial slope is Tr/v''^\8^). The curve may exist only if €i€3>l, which in 
the frequency domain corresponds to the cob<io:<^o:c/-\/2<^cop range. 

The above discussed results coincide, in the Xg/Xo <^<Cl region, with those given by the 
quasistatic approximation. 
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5. Brillouin Diagrams 



Keeping e^ and €3 as constant parameters in the Maxwell equations, it has been possible 
to solve the dispersion equation using a lengthy but straightforward procedure. Brillouin 

diagrams ( /5=y~ versus k = 2Tr/\o=o)/c, c being the light velocity j are, instead, more difficult to 

derive, because ei and €3 are functions of frequency. 

In order to obtain these diagrams from the previous section curves we must solve the 
following system: 



61H H 



^3" 






cog \ 



(15) 



For this let us keep constant the ratio «=— ? and choose all tlie pairs of ei and €3 values, 
wliich satisfy the relationship: 



ei = 



_ a(l — €3) — €3 
~a(l-e3)-r 



(16) 



Then tiie first equation of the system (15) provides, in the Brillouin kd versus ffd plane, a set 
of curves, €3, and a being constant for each curve. These curves will intercept the horizontal 
hues: 

M-M/Vl^3 (17) 

which represent the third equation of the above system (15), kpd being the plasma density 
parameter copd/c. For given values of the frequency independent parameters a and kpd, these 
intercepts at different 63 values describe the Brillouin curves: 



^d.=F{kd,a,kpd). 



(18) 



Typical Brillouin (Hagrams are plotted in figures 7 and 8 for kpd=l^.5 and for two a values, 
the first chosen in the a>l region (a=2.25), the second in the a<l region (a = 0.5). 

It is interesting to compare these curves with the corresponding curves obtained by 
Gould and Trivelpiece [1959] and by Snuillin and Chorney [1958] using the so-called quasi- 
static approximation. These authors have studied the more general case of a circular wave- 
guide partially filled with plasma; from their results the curves to be used in our case are easily 
derived by setting the waveguide diameter equal to infinity. 




Figure 7. Brillouin diagrarns for a=2.2r) and kpd = 
■3.5. 




Figure 8. Brillouin diagrams for a = 0.5 and 
kpd-'3.5. 
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No significant difference is found in the behavior of the Brillouin diagrams in the co<Cco^ 
region {kd<ikpd). In this region all the kd versus ^d curves start from the origin and approach 
asymptotically the value kpd or k^d, whichever is smaller, except the first solution mode 
which, when o^py-oob, approaches the asymptotic value kcd/'y/2={kld^-\-kld^)i/'\^. 

A behavior slightly more complicated than predicted from the quasi-static approximation 
is found in the co>cop region. Here the curves start from different points of the kd=l3d line 
between kpd and kcd and approach the asymptotic value kpd or k^d, whichever is larger. The 
slopes of these curves may be positive, negative, or partially positive and partially negative, 
so that the waves can be forward or backward. It is worth noting that in the quasi-static 
approximation all these waves are always backward. 

The following details can be added to the discussed features in the w^co^ region. 

(a) The starting points on the kd=l3d line are given from the simple relation: 



kd=kr>d 



V^+a+^ 



where 



p=[kpd/2X^^^\ 
(b) When co&^cop the wave of a given mode is always forward if: 

kld'-kld''>2(Xi''^y. 



(19) 



(20) 



If this condition is not satisfied, the wave may be backward, or start at low 13d values as for- 
ward and change to backward at large ^d values. 

For a plasma, which, instead of being stationary, has a uniform drift motion Uq along 
the 2 axis, the Brillouin diagrams are obtained from the above ones by substituting the Doppler 
shifted frequency co—I3uq to the signal frequency co. Our results are thus translated into the 
phase characteristics of space charge and transverse field waves of an electron beam. In 
this way they are useful for studies of traveling-wave and parametric amplification tubes. 

6. Power Flow 

The power flowing inside the plasma and the power flowing outside are obtained by 
integration of the Poynting vector over the respective areas. 
The z component of this vector inside the plasma is 

gP^^ (2tRV f N^Gl JUx^p/R) N2GI JKx2P/R) N,G,G, J,{x,p/R) Mx2p/R) 1 ^^i) 
where, for symmetry and simplicity, we have introduced the new symbols: 



^ L^l e/o(Xi) Xq i^o(^o) J 
V.X2 Jo{X2) Xo Ko{Xo)J 



^3 ^3^2 (Jl,2 



^Vi,2^ — ^1.2i^ — ■ ■ ^1 

€1 €1 a 



63(0:^— 1)6?,2— el 



N^=2gig2 



rr 10 ^3 €362 ffi + ^2 €3(ei— €3) g^ + l 

€1 €1 a ei a 



(22) 



The z component of the Poynting vector outside the plasma is 



l{x,) 



(23) 
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When integrations are performed, the ratio between the inside and outside powers can be 
written as: 



where 






Pa 



=(a)+(a+(a 



(24) 



(25) 



r Stpdp 

J R 

(Pv\ _N^ ( _ 2 JKxO I J\(x,) \ 
\Pa/i SWiMV Xi Joix,)^Jl{x,)J 

/PA _N^ /-. _ 2_ J^ , Jife)\ 
\PaJ2 SlgmV x,J,{x,yJl{x,)) 

(P,\^ 2N,{G,-G2) 
\PJ, {dl-SDg.g^M 

-[(My+<«-*''](-lit!-ill)-J 

The three terms of (24) are respectively proportional to the power transported by the 
field components of amplitude Ai, by the field components of amplitude A2 and by their mutual 
interaction (with reference to the expressions (5)). 

The Pp/Pa ratio is always real, also in the region where 8] and ^2 are complex; in fact, in 
this region Xi and X2, Qi and ^2, M ^^nd N2, Gi and G2 are all complex conjugates, so that the 
resulting {Pp/Pa)i and (PplPa)2 are complex conjugates too, and their sum is real, as it is the 
third term (PJPa)^- 




Figure 9. Power flow ratios Pp/Pa versus d/Xo in 
region I. 



1 I 


1 I -- r 1 

€,=-2 


T- 1 1 

€3=0.5 


~^ 


: SOL . \ 

A A 


n SOL. 


. \ L_. ., 




- 



2 4 6 ! 



16 1.8 2 



Figure 10. Power flow ratios Pp/Pa versus d/Xo in 
region III. 
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Figure 11. Power flow ratios Pp/Pa versus d/Xc in 
region IV. 
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The ratios P^jPa versus (i/Xo have been computed and plotted in fig*ures 9, 10, and 11 for 
the same ci and €3 values used in the numerical solutions of the characteristic equation, with the 
exception of the first solution in the €3<0 region. 

The ratios for this first solution have been evaluated in the limiting cases of zero and 
infinite rf/Xo values and have resulted much smaller than unity. The practical interest of this 
mode seems then very scarce; for this reason the complete curves have not been evaluated. 

A few more remarks seem worth mentioning: 

1. The Poynting vector has also a nonzero component, so that the power is actually 
flowing helically around the z direction. 

2. In region III the power flow inside the plasma may be backward. In this case, if 
Pp^Pa we have an '^average" backward wave, if Pp<^Pa an ^^average'' forward wave. The 
points labeled A in figure 5, where the group velocity is zero, correspond to the points A in 
figure 10 where Pp=Pa, so that, although power is actually flowing separatel}^ inside and 
outside the plasma, the total flow is zero. 

3. Also in region V the power flow in the plasma may be backward, although the total 
power is flowing forward, being PpPa- 

7. Attenuation Coefficient 

The attenuation coefficient a for the wave propagation is related to the power dissipated 
by losses per unit length P^ by the well-known formula: 



a= 



1 P. 



2 P -fP 



(26) 



We neglect the losses in the dielectric surrounding the plasma (usuall^^ air) and assume that 
the fields are those we have computed in the previous paragraphs, neglecting losses. In this 
case the losses are only due to the imaginary terms of the plasma dielectric constant and Pa 
is given by: 

Prf^Reair/co r E'D''2irpdp • (27) 

By means of (1) and (2), and recalling that Ep and E^ are 90° out of phase, this becomes: 

P,=o^eo r [e[\E,E^,+E^El)-j4\E,E;-E,Et)+e's'E,Ei]2Tpdp' (28) 

Jo 

Substituting the field expressions and performing the integrations, we finally obtain: 



where: 



'^-<■[d+')M"-^f"]-^'■£!-^■] 

'*'-■ [(1+0 s+l"-^i? »■>=<= SI-;"-] 






g\l'^Jl{x^)y9\\:^Jl{x.)] 9^92 4-4 ' J 



(30) 
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8. Conclusion 

Beside giving procedure^ foi' tlu^ solution of tlie characteristic equation of a uniform 
luagiieto-ionic plasma column of circular cross section and computing the X^/Xq and PplPa 
values for a few typical cases, our analysis has provided the following results, not given by 
the previous quasi-static analyses: 

(1) Propagation nui}^ take place in the cop<co<^coft region, provided (i/Xo is larger than a 
minimum cutoff value; in this case, at c?/Xo values sufficiently larger than the cutoff value, 
power flow takes place mainly in the plasma (Pp/Pa^l). 

2. In the co>aj6,co;, region we may have not only backward waves, but also forward waves. 
Backward waves carry most of the power inside the plasma. Here too propagation is possible 
only over a finite range of d/\o values. 

3. When co<Cco^(e3<^0), in the ranges where computations have been performed, power 
flow in the first circularly symmetrical mode takes place mainly outside the plasma, whereas 
all other circularly synmietrical modes show large Pp/Pa ratios. 



The authors thank A. Lemma, who carried out the numerical computations, applying 
also electronic computing techniques. The valuable lielp of Dr. A. Boschi is also acknowledged. 

9. Appendix A 

Regions I, IT, and III of figure 1 have in connnon the following features: 

(a) 0<63<1; 

(b) dl and 8l are always real; 

(c) ^2^ —1(52 versus a rises steadily from -co at a=l to —1 as a^oo); 

(d) — l:^r<0 (t versus a rises steadily from —1 at a = l to as a— >oo). 
Inspection of figure 2 shows that in this case: 

(a) When 6i>0 there are an infinite number of Xq solutions of (8); the mth solution falls 
between /"^^ (50 and /x'"^^ (5i) . As dl->+cx> this solution vanishes (xo^Xo"'^ /8i) and as dl->0 it be- 
comes infinite (xq'^XI"'^ /8i). In the r^Olhnit (a^oo) the Xq solutions become equal to i^^'"^(5i). 

(b) When 5i<0, no solutions of (8) exist. 

The above results determine the behavior of the Xq versus a solutions, for given €i and €3, 
if we know the characteristics of the 8^ versus a curves. These are as follows: 
(a) In region I, 8^ drops from +00 at a=l to at ao'. 



c^o=ei+V(e]-l)(ei-63) = ei + 62. (A-1) 

The parameter 8] becomes negative, when a>ao. 

(b) In region II, 8^ is always negative. 

(c) In region III, 8l is always positive and drops from +°° at a=l to 81, = — ^ >0 as a— >oo . 

10. Appendix B 

Regions IV and V, when: 

have in common the following features: 

(a) e3<0; 

(b) 5i and 8l are always real; 

(c) 5i and 8l versus a are respectively the upper and lower branches of a curve, the shape 
of which is that of a deformed parabola with the vertex at at and two horizontal asymptots 
for a->oo(8l^8'^ = -eJei and S'-^-l). 

(d) 0<r< 1 (r versus a drops from 1 at a^ to as a^oo). 
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Inspection of figure 3 shows that in this case: 

(a) When ^1^0, there is an infinite number of solutions. The smallest Xq (first solution) 
falls between zero and /^\8i); as 5i-^0, it remains then finite. All the other a^o solutions rise 
instead to infinity as dl-^Q, because the mth solution always falls between X['^~^^ /di and Xq'^V^i; 
if dl<iO. In the r=0 limit (a-^oo) the Xq solutions become equal to z^^'^^(6oo). In the t=1 
limit (a=ai, dl=8l), all the solutions go to infinity, except the first one which remains finite; 
in this limit (8) becomes an identity and the solution is given by the equation: 

[d{gF)/dSUa,^0. (B-2) 

(b) When — e3<5i<C0, there is always one solution alone. This solution falls between 
zero and v^^\di), and becomes equal to v^^\doo), as r-^O. 

(c) When 5i<^ — €3, there is one solution or no solutions depending on r being larger or 
smaller than: 

As T— >7o, Xq-^^ . 

Here too the Xq versus a behavior can be derived from the above results, whenever the 
characteristics of the Sj versus a curves are known. The basic features of these curves are 

(a) In region IV, 5L>0, so that there is always one region extending up to a=oo^ where 
5i>0. This region is the entire a>ai range, when: 

e,<e's=l(l-3e,-^l9el-10e, + l^ (B-4) 

and is the a>ao range [see (A-1)], when e^'^e^. 

(b) In region V, 51 <0 and then 5i<0 always. When 5t<~e3, no solutions exist as 
a-^00, because T=0<Cro', in this case, continuity and regularity suggest, as confirmed by 
extensive numerical computations, that no solutions exist at any a value. Propagation is 
then possible only if ei€3>l namely in the region: 



{ 



TSif"' (B-5) 

ei<e3<l/ei. 



1 1 . Appendix C 

Regions IV and V, when a^<ioi], have in common the following features: 

(a) 63<0; 

(b) 61 and 8] are complex conjugates; 

(c) r is a complex number, and \t\ =1. 

It is convenient to transform (8) into an equation of only real quantities. For this pur- 
pose, recalling that Bessel functions of complex conjugate arguments are complex conjugates, 
we write: 

8l2 = \8\' exp (±2je) 

gi,2=\9\ exp {±ji^) 

Jo(K2X^)=So{\8\xo,e) exp [±JM\^\xo,0)] 

Ji{8i,2X^)=S,{\8\xoe) exp [±JhW 3^0.0)]. J 
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(C-1) 



Equation (8) becomes: 

\8\' SVK'oixo) 



8\ SoKo(xo) L sin ^ "^ sin ?^ J* ^^ ^ 



When an Xq solution exists in the a^o^i region and it ends at ai with a finite value, then 
this solution has a smooth continuation in the a<C«i region. The Xq solution increases in tliis 
case as a decreases and becomes infinite at the a root of the equation: 

|gr+63 _ 63 cos (^+e)-cos (^-e) .^ ^. 

|5| ~ sin^ ^^ '^^ 
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